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Abstract. Multipartite entanglement, measured by the geometric entanglement (GE), is discussed for in- 
teger spin Valance-Bond-Solid (VBS) state respectively with periodic boundary condition(PBC) and open 
boundary condition(OBC) in this paper. The optimization in the definition of geometric entanglement can 
be reduced greatly by exploring the symmetry of VBS state, and then the fully separable state can be 
determined explicitly. Numerical evaluation for GE by the random simulation is also implemented in order 
to demonstrate the validity of the reductions. Our calculations show that GE is saturated by a finite value 
with the increment of particle number, that means that the total entanglement for VBS state would be 
divergent under the thermodynamic limit. Moreover it is found that the scaling behavior of GE with spin 
number s is fitted as a log(s-|- ^ -I-7) -1-5, in which the values of the parameters a, 13, 7, cr are only dependent 
on the parity of spin s. A comparison with entanglement entropy of VBS state is also made, in order to 
demonstrate the essential differences between multipartite and bipartite entanglement in this model. 



PACS. 75.10.Pq Spin chain models; 03.67.Mn Entanglement production, characterization and manipula- 
tion; 03.65.Ud Entanglement and quantum nonlocality 
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1 Introduction 

The crossover between quantum information and statistical mechanics fertilizes a distinct insight into the many-body 
effects and intrigues extensive interests [Ull]. Especially inspired by the pioneer works by Osterloh, et.al. and 
Osborne and Nielsen [4] , quantum entanglement in many-body systems has received great attention, especially on the 
critical systems [US]. This interest can be attributed to the following points. It is well known that the general properties 
of critical systems can be obtained by identifying the scaling behavior of the singularity for correlation functions, e.g. 
the correlation lengths and related critical exponents. Furthermore the vanishing of energy gap between the ground 
state and the first exciting state may lead to the long-range correlation even for infinitely separated particles at zero 
temperature. Thus it is conjectured that quantum entanglement, a special quantum correlation, could also display 
likewise behaviors at critical points. Some important progresses have been made in this direction. The recent studies 
for the block entanglement entropy have displayed the logarithmical divergency with the block size at critical points, 
and the universal area law has also been constructed exactly in one-dimensional spin-chain systems (see Ref. [5] and 
references therein). Furthermore the multipartite entanglement in many-body systems has also been shown to be 
maximal or to display sharp changes closed to the critical points P-(5!. 

From the point of quantum information, the long-range correlation implies that a stable quantum channel could be 
constructed easily for the transfer of quantum information |10| . For half-odd-integer antiferromagnetic Heisenberg spin 
systems, Lieb-Schultz-Mattis(LSM) theorem fT^ states that there is no energy gap between the ground state and the 
first excited state under thermodynamical limit, and the long-range order exists for any finite amount of anisotropy 
couplings. Then the long-distance quantum teleportation can be implemented in this system |101lll| . However the 
situation is different for integer spin-chain systems. Haldane's conjecture claims that there would be a finite excitation 
gap above the ground state for integer-spin antiferromagnetic Heisenberg chains, and then two-point correlation decays 
exponentially with the distance between two particles [13J. In order to verify this interesting issue, AfHeck-Kennedy- 
Lieb-Tasaki(AKLT) model was constructed, whose ground state is the so called valence-bond-solid(VBS) state |14j . 
And the two-point correlation functions for VBS state has been shown to decay exponentially because of the advent of 
the energy gap above the ground state [TS]. Consequently it is conjectured that quantum entanglement would be absent 
between two long distance separated particles in AKLT model. But the fact beats all. Verstraete, et.al. , demonstrate 
that the maximal entanglement between two infinitely separated particles can be founded even in the AKLT model 
by only imposing local operations on the rest of the particles [TB]. This finding enforces the reconsideration of the 
connection between quantum entanglement and correlation in many-body systems [17j[T8] . 

Although the absence of long-range correlation in AKLT model, den Nijs and Rommelse demonstrate that a 
hidden antiferromagnetic order still exists, defined by a nonlocal string order parameter(SOP) |19j . This nonvanishing 
SOP was then shown be related to the breakdown of hidden Z2 x Z2 topological symmetry in AKLT model |20) . 
Furthermore Venuti verifies the connection between the localized entanglement and SOP for spin-1 VBS state |18) . 
For higher integer-spin VBS state, several generalizations of SOP have been proposed to characterize the hidden 
topological symmetry pTH^ . 

Since it is well accepted that the local operation and classical communication(LOCC) cannot increase entanglement 
[24], the appearance of maximal localized entanglement in AKLT model thus means that quantum entanglement 
is predominant in VBS state. With respect of this point, the bipartite entanglement for VBS state has been first 
discussed [25,26 , where the block entanglement entropy is found to be saturated as 21n(s -|- 1) for spin-s VBS state. 
Unfortunately the behavior of entanglement entropy cannot present a clear explanation for the appearance of maximal 
localized entanglement in VBS state. Moreover the block entanglement entropy is insensible to the degeneracy induced 
by the topological symmetry. Furthermore the dependency on the parity of spin number s is absent for the entanglement 
entropy of VBS state, which has been manifested clearly by SOP [22^. In our point this defect could attribute to the 
fact that the overall information, embedded in many-body systems, is inevitably lost due to the trace-out of the 
superfiuous degrees of freedom for obtaining the reduced density matrix. And hence bipartite entanglement seems 
have limited ability of the complete characterization of VBS state. 

With these points the multipartite entanglement of VBS state is discussed in this paper in order to obtain the 
overall information embedded in VBS state. As shown in this paper, some distinguished features for VBS state can be 
obtained from the evaluation of multipartite entanglement, and moreover the multipartite entanglement of VBS state 
demonstrates different asymptotic behaviors from entanglement entropy. The paper is organized as the following. In 
Sec. 2, geometric entanglement(GE) is introduced for the measure of multipartite entanglement. And Valence-Bond- 
Solid state is defined with different boundary conditions in Sec. 3. In Sec. 4, the main part of this paper, we present a 
detailed evaluation of GE. Finally the conclusions and further discussion are presented in Sec. 5. 

2 Geometric Entanglement 

Geometric entanglement(GE) is used to measure the multipartite entanglement in this paper, which is first introduced 
by Shimony for bipartite pure state |27] and generalized to the multipartite case by Carteret et.al. [28', Barnum and 
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Linden ^5], Wei and Goldbart [30], to the mixed state by Cao and Wang [3T]. Geometric entanglement is a genuine 
multipartite entanglement measurement. The main idea of GE is to minimize the distance D between the entangled 
state and the fully separable state |^) in the Hilbert space, 

^= minim*') -|^>)|p}. (1) 

{|*>} 

Given the normalized \^) and |^), the evaluation of D is reduced to find the maximal overlap [50] 

Am ^ (2) 

Geometrically ^d??')) depicts the overlap angle between the vectors and |^) in Hilbert space. Then the larger 
^d*')) is, the shorter is the distance and the less entangled is In order to measure the entanglement in many-body 
systems, in this paper we adopt the definition in Ref. ^ 

e = lim , (3) 

L^oo Li 

where L denotes the chain length. Eq.Q actually defines the average entanglement per particle, and is also entangle- 
ment monotone. 

With this definition, the calculation for A(}^)) is generally difficult since the optimization. Recently two works 
prove independently that for the permutationally invariant entangled pure state the maximization can be attained 
necessarily by choosing a permutationally invariant pure fully separable state at the same time, whose amplitudes are 
all non- negative in a computational basis [32) . Moreover this conclusion has been generalized to any symmetric pure 
multipartite entangled state as shown in Ref. [33]. This point could be understood properly by noting the fact that 
K^l*')! — 1 means that and depict the same state, and both of them have the same physical features. While, 
— means that and have distinct physical features, and one can easily differentiate one from another 
by physical measurements. Thus for the purpose of maximizing it is necessary to find the fully separated |^) 

which has the same overall features to that of \^). 

Here we present a formal proof for this intuitive speculation. The start point for this proof is Cachy-Schwartz(CS) 
inequality 

\{v\w)\^ <{v\v){w\w), (4) 

where the equality occurs if and only if the two vectors \v) and \w) in Hilbert state are linearly related, i.e. \v) — c\w) 
for some scalar c. It should emphasize that c is unnecessary a constant, and generally would contain other variables. 
The trivial case for the equality occurring is that \v) and \w) are the same vector aside from the normalization 
constants. But for more general case, it is a nontrivial task for given \v) to find a physically different state \w) in order 
to maximize the overlap. However the task would become apparent if one rewrites the condition for the equality as 

\v){v\ \w) — {v\v)\w). (5) 

This transformation implies strongly that this task is reduced to find the eigenvector for operator This task 

seems trivial in the Hilbert space specified by vector \v). However it becomes nontrivial for enlarged Hilbert space. 
Fortunately recall that if two operators are commutative, they then share the same set of the eigenvectors, and then 
Iw) has same global symmetry to \v). Consequently this task can be reduced to find the all symmetries T for operator 
With this tricky, the optimal process in the definition of GE can be simplified greatly. 

3 Valence-Bond-Solid state 

Dependent on the boundary conditions, integer spin-s VBS state can be defined respectively as [22] 

L 

|VBS)PBC = n(4&l+i - 44+i)lO) (6) 
fe=i 

with periodic boundary condition(PBC) a^^]^^ — a'"i \ b^l]^^ ~ b^i \ in which a^^\ b'lj^ are the Schwinger boson operators, 
L is the chain length and |0) denotes the vacuum state, or 



L-l 



|VBS;p,g) = Qieft(al,6l;p) H (4^I+i " ^I4+i)' 



)Hght(aI,6[;g)|0), (7) 
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with open boundary condition(OBC) 

QHght(ai, 6[; g) = v/:^(ai)^-i(6[)^-«+i. (8) 

in which p,q — 1, ■ ■ ■ , s + 1 and „Cm denotes the binomial function. For OBC the VBS state is (s-l- l)'^-fold degenerate. 
VBS state can be rewritten as the matrix-product form [32] 

|VBS)pBc = Tr[5i • 52 • • • ffi] 

|VBS;p, q) = [gstart • .92 • • ■ 9L]{p,q) (9) 

where (p, q) denote the coordinates of matrix elements, and the matrix gi is (s + 1) x (s + 1) dimension, of which the 
element reads 

g,{p,q) = (-1)«+P- V.Cp-i,C,_i v/(s -p + q)l{s +p- q)\\s; q - p), 
5start(p, q) = ^/sCp-lsCq-l^/{s - p + q)l{s + p - q)l\s; q-p) (10) 

Obviously gi is related only to the i-th spin. Manifestly VBS state with PBC is pcrmutationally invariant, while it is 
not so for OBC because of the open ends. 

-Symmetry- VBS state is the ground state of Hamiltonian 

L 

i=l 

where the operator P^^i denotes to project the spins at sites i and i -\- 1 into the subspace with the total spin 

J = s + 1, • • • , 2s p^lTS] . Since Hamihonian Eq.(ITTl) is SU(2) invariant, then [H, n,^=i exp(i7rS'f )] = 0. This symmetry 
is also required for VBS state 

L 

[|VBS)(VBS|,^exp(^^^f)] = 0, (12) 

1=1 

where |VBS) denotes the VBS state without specifying the boundary conditions. 

-Magnetization- For PBC, each Sf eigenvalue occurs with equal probability |22| . and then the reduced density 
matrix of single spin is the unit matrix. This implies that each eigenvalue of the total 5*2 — X^i^i ^l^o occurs with 
equal probability, and the total magnetization {Sz)vbs in this case is zero. Furthermore {Sy)vBS — {Sx)vbs = 0. 
Thus VBS state for PBC is rotationally invariant. 

The situation becomes different for OBC. VBS is not necessarily rotationally invariant because of the open ends. 
Moreover the open ends may induce a local perturbation to the spin chain, and the hidden topological property for 
VBS state can be obtained from the system's responds to this local perturbation. 



4 GE in VBS state 

Given these properties, one can now determine the fully separable state \<P) = '^iLi \(t>i): in which \(f)i) — J2m=-s 1"^)* 
on the basis of {S'l'''; \'m)i\ and the coefhcient Cm is generally complex and dependent on the position i. With the 
requirement of Eq. fT^ . one has [|^)(^|,ni^i exp(z7rS'l*'')] = 0. Then two different situations can be identified as 

exp(»7r5«)|(/.,)f = \<\>,Y 

exp(z^^W)|^,)" = -|0,)" (13) 

where 

[s/2] 

E 

fe=-[s/2] 

[(s-l)/2] 

E 4'^|2fc+l)» (14) 

fe=-[(s+l)/2] 
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Fig. 1. (Color Online)The randomly numerical simulation of — K'^|vbs)| VBS state with PRC and spin number s = 1, 2, 3 
respectively. There is 1000 sample data for every simulation. The plots of (a), (b) and (c) correspond to the simulations without 
any other restriction on the coefficients c^'' of the fully separable state except of the requirement of normalization. While 
for plots (d), (e) and (f) the requirement of permutational invariance for |#) has been implemented, but Cfe for different k is still 
chosen randomly. The (blue) solid line corresponds the value of GE after applying the all approximations stated in the paper. 



and [n/2] denotes the maximal integer number not bigger than n/2. Then one has to find the maximal overlap among 
{|(VBS| (H)f^i |(VBS| (Kif^i Consequently GE is determined by the minimization 

(T = min{ep, £„} (15) 

where ep(„) denotes GE evaluated with the separable state <8)fili|0)f^"^ respectively. Dependent on the boundary 
conditions, the following discussion has to be divided into two subsections. 



4.1 Periodic Boundary Condition 



As displayed in Sec. 3, VBS state is permutationally invariant in this case. It has been proved in Refs. that 
the fully separable state is necessarily also permutationally invariant in order to maximize the overlap. Then the fully 
separable state \<1>) becomes in this case 



(16) 



Another important feature is that VBS state is rotationally invariant for PBC, and hence has zero total spin |15j . 
It is natural to endow the fully separable state |^)p(n) with the same character. This point is evident to note that it 
belongs to completely different space for the states which has vanishing total spin or not, and then the overlap between 
states in the two distinct spaces must be zero. 



With this requirement, the coefficient c^.*"* preferably has the same amplitude, independent both of the position i and 

quantum number k. Moreover c^^^ should preferably be in phase in order to maximize the overlap with respect of the 

gi has a block-diagonal 
of 



special form of VBS state shown by Eq.(IS]). This point becomes clear by noting first that ^ 
form dependent on the parity of p, q for any spin s (see Appendix A for more detail) . Set Ck — o. 
depend only on the difference |g — p|. And then 



which Oq^p and 



v-g 



^"-pgp^qleven p and q)\ 



p(7p_<j(odd p and q)\ 



(17) 
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where the dashed boxes denote the sub-matrixes satisfying the conditions stated in the brackets, and the element of 
the sub-matrix is read 



Then one has 



\pmVBS)v^c\ = |Tr[p(<Z'|gi • 52 • • • = 

s+l 



s+l 



pi:P2,---pt;9i,<32.--- ,gi, = i j=i i=i 
L 



n 



pi-P2,---PL;qi,q2-/-- ■qL=i i=i 

pi,P2,---PL;qi,q2,--- ,qL=i i=i 
s+l L 

pi,P2,---PL;qi,q2,--- ,qL=i i=i 



where the following relation is used 



ni 



< 



Pi-qi 



(18) 



(19) 



in which the equality occurs if and only if \cp.-q- \ = \cp-^q. | = |cp| for arbitrary i 7^ j. 
With even L and noting the common factor (— 1)*+p^i in gp q, then 



s+l 



I 



L 

pi,P2,---PL;qi,q2,--- ,qL=i i=i 



(20) 



which means that one need only choose Cp being non-negative real for the maximization of the overlap. As for odd L, 
one should note that gi is traceless and the overlap always vanishes. 

For 10)", "(01 (E) gi has a block-anti-diagonal form with opposite sign for the two blocks, as shown in Appendix A, 
and thus the overlap always vanishes for odd L. For even L one only notes that the product of two block-anti-diagonal 
matrixes is block-diagonal, and the two diagonal block has the same sign after this product. Then the analysis is 
applicable in this case. 

Then the fully separable state can be determined exactly 

L ls/2] 

i=l k=-[s/2] 
L [(^-l)/2] 

|<?)„-(g)c„ J2 |2^+1)» 

1=1 A;=-[(s+l)/2] 



(21) 



in which Cp — l/-\/l + 2[s/2], c„ ~ l/-\/2[(s + l)/2]. Then with respect of the matrix- product form of VBS state, the 
overlap can be expressed readily as 



An = 



\pmVBS)pBc\ 



Tr[p(<^|.gi-g2---.9L] 



^PBc{VBS\VBS)pBC VT^y[9*i-9*2---9l<»9i-92---gL] 

\n{^\VBS)pBc\ ^ Tr[„(j>|.gi - , 92- --gL] 

^PBc{VBS\VBS)pBC VT^r[9l-9*2---9l<^9i-92---9Ly 



(22) 



where g* denotes to take Hcrmitian conjugate for the matrix elements, but not transpose the matrix. 

In order to demonstrate more explicitly the validity of the statements above, the randomly numerical simulation 
has been implemented for several simple cases, as shown in FigH] For FigUJa), (b) and (c), there is no any other 
restriction on the choice of c[.''' except of the requirement of normalization. It is obvious that the sample data are 
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Fig. 2. (Color Online)(a) versus the particle number L with PBC. In this figure the Arab numbers denote the spin and 

p,n denote the plotting for ep,e„ respectively; (b) the scaling of Sp with respect to spin s and the fitting function /(s); Use 
/(s) to fit GE for even(c)and odd(d) spin-s case. 



Table 1. The chosen parameters for the fittings of GE in Figsl^c) and (d). 





a 


/3 


7 


5 


Fig. Etc) 


1.41 


-0.39 


2.82 


0.83 


Fig. [ad) 


1.15 


1.67 


-1.83 


1.33 



always bigger than the value of GE evaluated after applying the approximations. Furthermore the approximation 
that is a positive constant independent both of the coordinate i and quantum number fc, has also been verified 
numerically as shown in Fig|ljd), (e) and (f). The results show that this approximation is really proper and correct. 
With the assistance of the numerical simulation, GE can be evaluated readily. 

For arbitrary spin-s the analytical expressions for the overlap are difficult to obtain. Thus ep,en are plotted 
numerically in Fig. [D^a). It is obvious that ep,en tend to be saturated with the increment of particle number L. 
Moreover one can note Sp < e„ under large L. As our study is implemented mainly under the limit of large L, the 
discussions below would focus on Sp. 

It should emphasize that since £p(„) denotes the average GE per particle, the saturation for ep(„) shown in Fig. 
[2ta) means that the globe GE Eg for VBS state, defined first in Ref. [7], should be proportional to the particle number 
L, and is divergent under L — >■ oo, i.e. 

Eg = SpL. (23) 

This phenomena is distinct from the entanglement entropy of VBS state, which is shown to be saturated as 2 ln(s + 1) 
because of the appearance of energy gap above the ground state !25''26 . However the divergency of Eg strongly implies 
that there is hidden degeneracy not captured by entanglement entropy. With the findings of the breakdown of hidden 
Zs+i X Zg+i topological symmetry [lOlllS], it hints that GE would be also sensible to the existence of topological 
phases in many-body systems. A further discussion would be presented in the final section. 

Furthermore, as shown in Fig|5][b) where the scaling behavior of Sp with spin-s is plotted, £p shows two different 
scaling behavior, dependent on even or odd spin s. This feature is consistent with the observation that one have to 
define different SOPs for odd or even s VBS state. A fitting function /(s) can be constructed for identifying the scaling 
behavior of GE with spin s, 

/(s) = alog(s-H--t-7)+(5 (24) 
s 

in which a, /3, 7 and a are tunable parameters. The chosen values of these parameters for the fittings of Sp are listed 
in Table HTTl It has been shown that the scaling behavior with spin s for entanglement entropy of VBS state pSlf^ is 
attributed to the tendency to be unit form for the reduced block density matrix under limit L — >■ 00 j26j . Comparably 
it is unclear until now for the distinct feature of GE disclosed by Eq. ([M)) . 

The differences between entanglement entropy and GE in VBS state can intuitively attribute to the trace-out of the 
superfiuous degrees of freedom for obtaining the reduced density matrix in order to calculate the block entanglement 
entropy. Thus the overall information about VBS state is lost because of this trace-out. With respect of this point, it 
seems not enough by only sampling a portion of the system to obtain the overall information of many-body systems. 

Another interesting character for this numerical evaluation is that Sp^En are sensitive to the parity of particle 
number L. As for odd spin s case, Ap is always vanishing for odd L and then Sp is infinity. This point can be easily 
verified (see Appendix A). Whereas for even spin-s our calculation shows that Ap for odd or even L tends to be 
consistent numerically with the increment of particle number. The similar behavior also happens for the evaluation of 
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s=1, E =1.58496, L=1000 




Fig. 3. (Color Online)Color Online)Tlie randomly numerical simulation of — ^ {■^[v'bs) | VBS state with OBC and spin 
number s = 1, 2, 3 respectively. There is 1000 sample data for every simulation. The plots of (a), (b) and (c) correspond to the 
simulations without any other restriction on the coefficients c^''' of the fully separable state \<P) except of the requirement of 
normalization. While for plots (d), (e) and (f) the requirement of permutational invariance for has been implemented except 
of the end pstart, but Ck for different k, including that for pstart, is still chosen randomly. The (blue) solid line corresponds the 
value of GE after applying the all approximations stated in the paper. 



4.2 Open Boundary Condition 

The situation becomes complex for OBC. Although the requirement Eq. (jl2l) persists, the permutational invariance for 
VBS is broken in this case. Thus the assumptions made in the preceding subsection has to be reconsidered seriously. 

The crucial step is still to determine the overlap |(VBS;p, It is interesting to note from the matrix-product 
form of |VBS;p, g) that except of the end ijstart the remaining giS of this expression are obviously permutational 
invariant. With respect of the findings in Ref. [31], one could reasonably assume 

{VBS;p,q\'P) = [start ((/''|.gstartl(0|.gi ■•• L-l(0|.gL-l](p,g) (25) 

Now it is crucial to find the relation between |0') and |(/)). 

As stated in Ref. [35] , the boundary effect vanishes exponentially with increasing length of the system if one focuses 
on the local operator. From the definition of GE Eq.([2]), it actually imposes a locally independent operator \(j)i){ipi\ on 
the individual spin since the choice of is independent with each other in the definition. Thus it is an acceptable 
approximation to ignore the effect of f/start on determining |0') under the limit of large L, and let |0') — {(j)). Another 
important character for ^start ■ 92' ■ ■ dL is that the eigenvectors of single spin 5'| occur with equal probability |22) . 
which also means that the different boundary condition represented by (p, q), happens with the same probability. With 

these observations, it also is a reasonable assumption that the coefficient cj. could be same to that for PBC. Hence 
for OBC the fully separable state is same as that for PBC (see Eq. (HI])). 

In order to demonstrate the validity of these approximations above, the numerical simulations are implemented 
by randomly choosing c'j^\ Similar to the simulation for PBC, the plots FiglHJa), (b) and (c) display the results of 

simulation without any other restriction on c"^^ except of the requirement of normalization. For plots FiglJl^d), (e) 
and (f), the permutational invariance is imposed on the other particles except of the end (jstart, in order to verify the 
approximation that c[.'' is a positive constant independent of the coordinate i and quantum number k. These plots 
manifestly display the validity of this assumption. Furthermore the effect of OBC has been considered independently, 
as shown numerically in Fig|3J For this simulation, the approximations stated in the paragraph above are implemented 
for the particles except of the end, in order to highlight the effect of boundary condition. Obviously with the increment 
of particle number L, the effect of OBC tends to be vanishing. And so except for s = 1, this observation means that 
these approximations are valid even for very small L. As for s = 1, our further numerical simulation shows that with 
the increment of L, all data would become consistent with the result obtained by the approximations above. It should 
point out that with the consideration of the small L e = e„ , which is obviously below Sp as shown in Fig|5l For the 
convenience of the comparison with previous observations, only Sp is plotted in Fig|31[al) and (a2). 
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Fig. 4. (Color Online) The numerical simulations for the effect of OBC with different L. Except for the end (?start, the state 
\(j>} is chosen based on the approximations stated in the paper, and |(;i')start is still randomly chosen. There is 1000 sample data 
for every simulation. 



With these approximations, it is not difficult to find the overlap |(VBS;p, q\'P)\- In order to eliminate the dependence 
on special boundary condition, the average for all degenerate states is adopted 

l^(|VBS))| ~(7TIF^,(VBS;p,,|VBS;p,<z)' ^^'^ 

with 

L-l 

{VBS; p,q\<Pf = Lgstartl'/') [] aM)^]^- (27) 



Since the elements of gi and \(j))i are real, one obviously obtain under the limit of large L 

L-l 

|^(|VBS))p = _Tr(2)[g,t„t|0) J] 5,10),]. (28) 

where Tr^^-* [M] denotes the sum of the square of all the elements in matrix M, and the normalization constant 
CL^{s + for large L 

For s = 1 one has |yl(|VBS))p = for separate states |<?)p and |<?)n. Then we have Sp = e„ = logj 3 under 
L oo. In a recent work of Oriis [5], the author discussed the block GE for spin-1 VBS state. Our result is consistent 
with Oriis' result when the block has only one particle. For s = 2, |yl(|VBS))p = [(4/3)^(1 -fV6)^ + (l-\/6)^-|-4^]/cL 
for \'P)p, and 3^" /cl for separate state |<?)n- For higher s, the exact expressions for yl(|VBS)) is difficult to find. Hence 
our discussion below relies on the numerical evaluation. 

In Fig[5]ep,£„ are plotted respectively. As shown in FiglSI^a), both £p,£n tend to be saturated with the increment 
of L, similar to the case of PBC. Moreover Sp and e„ are also generally unequal, as show in FiglSJa). An intricate 
phenomena is that because of the finite-size effect Sp > En for small L, while under the limit of large L, Sp < En- 
Comparing Fig[S]and Fig|5J our calculation show that the two Sp are equivalent numerically under large L whether for 
PBC or OBC, so does for e„. While for small L, the value of Sp is obviously dependent on the boundary conditions. 
These findings are consistent with our assumption that the boundary effects have exponentially decayed affect with 
the increment of system size. It means that for VBS state, geometric entanglement under thermodynamic limit is 
independent on the boundary condition. Hence the scaling behavior of Ep with spin s shown in Fig. [5{b), is same to 
Figl2Ib) and has the same fitting functions. 
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Fig. 5. (Color Online) (a) £p, £„ versus the particle number L with OBC. In this figure the Arab numbers denote the spin and 
p, n denote the plotting for ep,£„ respectively, (b) the scaling of with the spin s, which is very similar to Fig[2jb). 



5 Conclusion and Further Discussion 



In conclusion the geometric entanglement is discussed explicitly for integer spin Valence-Bond-Solid state in this paper. 
In order to reduce the optimization in the definition Eq.Q, we point out that the overlap between the entangled 
state and the fully separable state \^) is dependent on the common symmetry shared by two states. With this 
important observation, the fully separable state which maximizes the overlap with VBS state, can determined 
explicitly whether for PEC or OBC. For PEC, |«?) is permutational invariant and can be determined exactly. While 
for OBC, since the perturbation induced by the boundary condition, the fully separate state |^) can be decided when 
the system size tends to be infinite, and the effect of boundary condition can be ignored in this case. Furthermore 
in order to display the validity of the approximations made to determine |^) exactly, the numerical simulations have 
been implemented, which demonstrates clearly the validity of the approximations. 

An intricate property for the evaluation of GE is that ep(„) is dependent on the parity of the particle number L. The 
similar feature can also be found for the evaluation of entanglement entropy in one-dimensional system |36j , however 
in which case the dependence of the parity of particle number is because of the boundary conditions. While for VBS 
state, the parity dependence of GE is mainly because that the elementary matrix (g) gi is traceless. Another 

interesting feature for GE is that the for PBC and OBC coincides under large L. Compared to the entanglement 
entropy, which heavily depends on the boundary conditions 151, this feature can be considered as a manifest of the 
overall character of multipartite entanglement, and shows that the local boundary condition would be negligible when 
one focuses on the global entanglement in the system. 

The important observation for GE is that becomes saturated with the increment of particle number L as 
shown in Figj^J^a) and FiglSJa). Since denotes the average entanglement per particle, this feature strongly implies 
that the total entanglement for the whole system is proportional to the particle number L, and the globe GE Eq is 
divergent under L — ^ 00. Given the appearance of the finite energy gap for the system described by the hamiltonian 
Eq. pT|) [T5] , this divergency means that there is hidden symmetry for the hamiltonian EqUT] Together with the finding 
of the hidden Zs+i x Z^+i topological symmetry for Ea. lfTTj) [IHIISS], it implies that GE could be used to detect the 
topological phase transition. Moreover this divergency also present a clear explanation for the appearance of maximal 
localized entanglement in AKLT model (TB] ■ Since the degeneracy of ground state of Eq. (|TT|) still exists manifested by 
the divergence of Eq under X — >■ 00, it means that the long-range correlation could still be founded between any two 
particles. Another important observation is that the scaling behavior of Ep with spin number s is obviously dependent 
on the parity of spin s shown in Figl^fb) and FiglSlJb). This feature is consistent with the observation that different 
SOPs have to be chosen for odd and even spin-s VBS state in order to identify the different topological symmetries [25]. 

Given these points it is believed that GE could include more information of VBS state than that of entanglement 
entropy. These critical differences between GE and entanglement entropy for VBS state suggest that multipartite 
entanglement would be more popular for our understanding of the many-body effects since GE presents an overall 
description for many-body systems. Moreover it also hints that GE could be used to mark the different topological 
phases, of which the identification is still a difficult task. Another intricate points is that the convergence of GE is 
much slower than that of entanglement entropy 37. ■ In our own point, this feature can attribute to the long-distance 
correlation characterized by SOP. 

Finally, it should point out that the definition and measurement of multipartite entanglement are a difficult task 
in general [38]. By the discussion of GE in VBS state, we should express our own point that it is efficient for the 
construction of a comprehensive understanding of multipartite entanglement to found its connection with the diverse 
physical phenomena, especially in many-body systems. Hence we wish this work may be useful for this purpose. 
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Appendix A 

This point can be verified by noting 



T^[p('?l5i-.92---gL]=Tr[n^(0l®3d 



(29) 



1=1 



With respect of the expressions for gi (Eq. ([TU])) and ^{(pl (Eq. (pTj) ). the relation above can be rewritten as a block- 
diagonal matrix for s = 2k — l(fc = 1, 2, . . . , n) 



' 9^ 



91,1 91,3 

93,1 93,3 



ffl,2fc-l 
93,2k-l 



52fe-l,l 52fe-l,3 ■ ■ • .g2fe-l,2fc-l 







92k,2k <72fe,2fe-2 ' ' ' 92k-2.2 
92k-2,2k 92k-2,2k-2 ' ' ' 92k-2.2 



g2.2k 



92,2k-2 



92.2 



(30) 



where 



9p,q - {-iy+^-^^sCp-isC,-i^{s~p + q)\{s+p- q)\. 



The first submatrix correspondes to the case both of p and q being odd number, while the second submarix corre- 
sponding to even p, q. Noting that 2k-iC2k-2m-i = 2k-iC2m+i-i, one conclude that the two submatrices are identical 
execpt of a negtive sign. Moreover this matrix is independent of the coordiante i of particles, and then for odd L the 
trace of Eq. ((29|) is always zero. 

However although the block-diagonal form still exists for even spin s, the two diagonal blocks do not have the same 
dimensions. The conclusion above is not correct for even s VBS state. 

The analysis is similar for the evaluation e„. In this case ® gi can be formulated as a block-antidiagonal matrix, 

i.e. 



(p\®gi^ Cn 









52,1 32,3 ■ ■ 


52,2j-l 1 


54,1 54,3 •• 


54,2j-l i 


52i,l 52i,3 ■ ■ 


52i,2j-l i 



(31) 



V 

where j = [(s + l)/2], j = [s/2] + 1 (the function [/] denotes the integer part not more than / ). Then for odd L 
'^[11^=1 ® 9i\ always is zero. 
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